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Previous studies have shown that small changes can be monitored in a scattering medium by observ- 
ing phase shifts in the coda. Passive monitoring of weak changes through ambient noise correlation 
has already been applied to seismology, acoustics and engineering. Usually, this is done under the 
assumption that a properly reconstructed Green function as well as stable background noise sources 
are necessary. In order to further develop this monitoring technique, a laboratory experiment was 
performed in the 2.5MHz range in a gel with scattering inclusions, comparing an active (pulse-echo) 
form of monitoring to a passive (correlation) one. Present results show that temperature changes in 
the medium can be observed even if the Green function (GF) of the medium is not reconstructed. 
Moreover, this article establishes that the GF reconstruction in the correlations is not a necessary 
condition: the only condition to monitoring with correlation (passive experiment) is the relative 
stability of the background noise structure. 



PACS numbers: 43.40.Ph, 43.60.Tj, 43.40.Le, 43.20.Fn 

I. INTRODUCTION 

In order to image a complex medium the impulse 
response, or Green function (GF), of that medium is 
needed. Classically, the impulse response is retrieved by 
active means: a signal generated at one point {eg. an 
earthquake) is recorded at another (a passive receiver 
hke a seismic station), and this record is treated as the 
band-pass filtered GF. Over the last fifteen years, devel- 
opments in helioseismology^ and in acoustics^ showed 
that information about a medium can be extracted 
from diffuse (coda) waves or ambient background noise. 
Since then, seismologists are turning to this passive 
imaging technique^. In this latter case the seismic coda^ 
or seismic noise^ is correlated to reconstruct the GF, 
either by averaging over space, time, or both. Passive 
imaging requires some assumptions: one needs uniformly 
distributed noise sources and/or long enough record 
duration to get the correlation function to converge to 
the GF. 

Monitoring dynamic media is a separate issue and, 
as we will see, is based on different (and weaker) 
assumptions. In the 80's, Poupinet et al.— proposed to 
use coda waves to monitor velocity changes in scattering 
media (the after- mentioned doublet technique). At first 
glance, coda waves might appear as a jumbled mess 
of wave arrivals. In fact, they consist of the waves 
which have described long, scattered paths through 
the medium, thereby sampling it thoroughly. As a 
result, these scattered waves are more sensitive to small 
variations than the ballistic waves. For this reason, the 
information in the coda can be exploited to monitor 
small changes in a medium. This technique, analogous 
to Diffuse Wave Spectroscopy (DWS^) in optics, was 
later named coda wave interferometry (CWI^). It tracks 



the tiny phase changes in the coda that are caused by 
velocity changes in the medium. A major issue of the 
doublet technique is that it requires stable reproducible 
sources, which are hardly available in seismology. Thus, 
a more recent idea was to combine noise-based passive 
imaging with the doublet tcchniqueiSiiii'i^iii^. First, one 
correlates the background noise between two receivers. 
Second, one analyzes small phase changes at large lapse 
times (coda) in the correlations. This forms the basis of 
passive monitoring (or passive image interferometryii) 
with seismic noise. Noise based passive monitoring 
seems to simultaneously require two conditions. First: 
a homogeneous distribution of sources in space, and 
second: temporal stability of these sources. 

In seismology, most of the noise (between 0.01 and 
1 Hz) generated in the ocean o^'*'^^ shows strong spatio- 
temporal variabilities. This feature is in favor of passive 
imaging as long as records are long enough (duration of 
the order of a year) to average over a large distribution of 
sources. But to passively monitor dynamic phenomena 
over a few days or less, this feature seems very unfavor- 
able. In this paper, we investigate the effect of these 
non-ideal conditions on the reliability of passive moni- 
toring. We will also examine if the GF reconstruction in 
the correlation is a necessary condition to perform pas- 
sive monitoring. To that end, we test the passive mon- 
itoring technique under degraded conditions in a con- 
trollable (laboratory) environment. In section II of the 
present paper, we describe the experimental setup and 
our motivations to do small-scale seismology. In section 
III, we compare two different data processing procedures 
to extract velocity variations from the coda wave. One 
is the doublet technique introduced twenty years ago. 
The other, referred to as stretching, was developed more 
recentlyiiii^iii. Advantages and drawbacks of both pro- 
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cedures are discussed. We also investigate the robustness 
of these procedures when noise is introduced in the sig- 
nal. In section IV, we investigate if passive monitoring is 
still possible when the GF is not properly reconstructed 
in the correlations. Finally, in section V, we test the ro- 
bustness of passive monitoring in the case of temporally 
changing distribution of sources. 

II. METHODOLOGY 

A. Motivations for doing analogous laboratory experiments 





Seismology Ultrasound 


Wavelength 


km 


mm 


Frequency 


mHz - Hz 


MHz 


Total size 


10^ km 


m 


Total duration 


month-year 
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TABLE I. Comparison of the physical parameters between 
seismology and ultrasound. 

Seismology is based on the observation and process- 
ing of natural vibrations. In a passive field experiment 
where seismic waves originate from earthquakes, scien- 
tists are facing two simultaneous problems. They neither 
know the source location with sufficient precision, the 
source mechanism nor the medium of propagation. It 
is therefore very complex to image the source and the 
medium at the same time. By reproducing some fea- 
tures of the seismic propagation in the lab and employ- 
ing controlled sources and sensors, we can focus our ef- 
forts on the physics of the wave propagation and develop 
new methods more comfortably. In laboratory-scale seis- 
mology, we control for instance the size of the medium, 
the scattering properties, and the absorption. We are 
then able to adjust one parameter at a time and test 
the physical models and imaging techniques we develop. 
But the main reason for carrying out analogous ultra- 
sonic experiments is more tactical: it is related to the 
order of magnitude of the physical parameters as recalled 
in Table HI Ultrasonic wavelengths are on the order of a 
millimeter, meaning that experiments are physically easy 
to handle. Additionally, the duration of a single ultra- 
sonic experiment is very short (one minute) compared to 
seismology where we have to wait for earthquakes (year) . 
This characteristic allows us to achieve many more exper- 
iments in the lab, and test many parameters over a wide 
range of magnitudes. In the view of testing processing 
technique to monitor weak changes, it is also of first im- 
portance to perfectly control the origin of the change in 
the medium. This is quite convenient in the laboratory, 
but almost impossible in a natural environment. These 
are the reasons why several seismology laboratories have 
decided to develop analogous experiments for method- 
ological developments^'^*'^^'^". Our article presents one 
analogous ultrasonic experiment not only devoted to the 
study of the physics of wave propagation in heteroge- 
neous media, but also to the development of new tech- 
niques applicable to seismic waves in geosciences. 



B. Scattering properties of the medium 

We perform the experiment on a 80 mmx64 mm block 
of Agar-Agar gel which consists of 95% water and 5% 
Agar (by weight). 8.5% of the volume of the gel con- 
sists of small air bubbles, with diameters between 100 
lim and 1 mm. These bubbles render the medium mul- 
tiply scattering. The source emits a pulse at 2.5 MHz 
(100% frequency bandwidth). For simplicity, we neglect 
the electronic noise in the experiment. Since shear waves 
are strongly attenuated, we assume that only P-waves are 
propagating in the medium and are eventually recorded. 
To estimate the scattering properties of the medium, we 
performed several experiments in the transmission config- 
uration for several medium thickness. From the atten- 
uation of coherent plane wave, we obtain an estimation 
of the elastic scattering mean free path averaged in the 
working frequency band: £e ~ 3.5mm. Since the scat- 
terers' size is smaller or of the order of the wavelength, 
scattering is isotropic and we expect a transport mean 
free path £* of the same order. In Fig. [T] we plot an 
example of a diffuse record transmitted through 64 mm 
of our heterogeneous medium. A theoretical fit is ob- 
tained from the two-dimensional (2D) diffusion equation 
(including reflections from the sides) and plotted as a 
broken black line. The diffusion constant is D = vp£* /2 
and we assume ^ £*■ The absorption length £a is the 
fit parameter, the best fit is obtained for £a = 200 mm. 
This corresponds to an absorption three times stronger 
than in pure water— which is due to dissipation by the 
agar material. 




250 500 750 1000 

Time (|i s) 



FIG. 1. Gray line : acoustic field transmitted through the 
64 mm thick bubble-gel mixture (in normalized amplitude). 
Black broken line: 2D diffusion equation for £* — 3.5 mm and 
£a = 200 mm. 



III. COMPARISON OF DATA PROCESSING 
TECHNIQUES 

Two processing techniques have been proposed in the 
literature to estimate relative velocity changes dV/V 
in the diffuse coda. The first one, called the seismic 
doublet technique, was developed for geophysical pur- 
poses about twenty years agoi. The idea is to measure 
a time-shift between two different records in limited 
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time-windows centered at t in the coda. By repeating 
this procedure at different times t, it is possibie to piot 
the delay St versus t. The velocity variation is simply 
the average slope of St{t): dV/V = —St/t. Doing so, we 
implicitly assume that the time-shift is constant within 
the considered time-window, which might be not the 
case. This processing found remarkable applications 
in geophysics, including recent developments in vol- 
cano eruption predictiortiS: and active fault monitoring^^. 

Another ide a^^'^^ is to interpolate the coda at times 
t{l — e) with various relative velocity changes e. This 
corresponds to stretching the time axis. The actual ve- 
locity change is obtained when the interpolated coda best 
fits the original data. Because we do not assume a con- 
stant time-shift in the considered time window [0 T], we 
can process the whole data at once, which is expected 
to result in a more stable, and thus more precise, estima- 
tion of dV/V. One drawback is that this latter processing 
assumes a linear behavior dt{t) versus t, or a constant rel- 
ative velocity change dV/V — e, which is sometimes not 
the case in complex heterogeneous media. No quantita- 
tive comparison between these two techniques have been 
established in the literature. In the following section, we 
propose to test both techniques on the same data set, 
and analyze their sensitivity to the SNR of the records. 

A. Active experiment: high quality data. 

In this experiment, we attach a set of transducers on 
one side of the gel which act as both sources and receivers 
of the signal (figure [2]-left). The source emits a 2.5 MHz 
pulse. The signal is collected on the same transducer 
(R) in the pulse-echo configuration. This procedure is 
repeated on 7 different channels. As the gel contains 
a large amount of scatterers, the emitted waveform is 
multiply scattered before reaching the transducer again. 
A typical ultrasonic record is plotted in Fig. [1] Note 
that the early 5fis are muted for technical reasons. This 
record is composed of the GF of the air-gel mix sample 
G{R,R,t) and the source wavelet e{t): 

ho{t) = Ga{R,R,t)®e{t) (1) 

where <Si stands for convolution. This experiment is re- 
peated 4 times while the temperature of the medium 
slowly increases by about 0.8°C, as measured by a dig- 
ital thermometer placed beneath the gel. We assume 
that the first effect of a temperature change is to stretch 
the record in time by £fc = dV/V, and to additionally 
slightly distort it^^. This weak distortion, noted f{t), is 
not studied here, although it contains precious informa- 
tion about the medium and its evolution. An example 
of two records is displayed in Fig. [31 After a (small) 
temperature change, the record rewrites: 

hk{t) = Gk{R,R,t)(g)e{t) (2) 
^[GoiR,R,t{l+ek)) + fit)](g>eit) (3) 



Active configuration Passive configuration 



— 64 mm - 




FIG. 2. Experimental setup. Left: the active experiment in 
the pulse-echo configuration. Right: the passive experiment 
in the transmission configuration before auto-correlation. 




10 20 30 40 50 

Time (n s) 



FIG. 3. Example of two records ho{t) and /i3(t) acquired 
at the same position (same transducer) but at two different 
dates. Between the two acquisitions, the temperature has 
increased by 0.8° C, which is hardly visible in the early part 
of the record (inset between 10 and 11 /is) but very clear in 
the late coda (inset between 42 and 44 /is). 



For each temperature fc, the record hk{R, R, t) is com- 
pared to the reference waveform hQ{R, R,t) to evaluate 
the relative velocity change in the gel sample. Two pro- 
cessing techniques have been proposed in the literature to 
estimate dV/V: the doublet technique and the stretching 
technique. 

1. Doublet technique 

The doublet technique also known as cross-spectral 
moving- window technique (CSMWT^^), computes the 
phase shift between records for consecutive, overlapping 
time windows. For a given window, the time shift is 
assumed to be constant and is estimated in the frequency 
domain by measuring the Fourier cross-spectrum phase. 
This estimator uses an accurate, unbiased Wiener filter 
technique^! and produces an estimate whose confidence 
interval is controlled by the coherence values in the 
frequency band used for the analysis. The method 
can then measure arbitrary time-shifts between two 
records with enough similarity (or coherence). The 
key parameter in this analysis is the Fourier transform 
window length. The length choice is a trade-off between 
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shift estimate accuracy, and the time resolution of 
possible temporal variations. 



We use this doublet technique to compute the time 
shift between the records. The time shift between the 
two different records is measured in the coda between 
12.5 and 50 fis. If the velocity changes homogeneously in 
the medium, the propagation time will vary proportion- 
ally to the propagation distance, producing a phase shift 
between records varying linearly with lapse time. The 
relative velocity change can be retrieved by measuring 
the slope of the phase shift as a function of lapse time, 
as shown in figure [51 



FIG. 4. Delay time evaluated for different lapse times in the 
coda from the doublet code. The broken line is the linear 
trend whose slope yields dV/V between fe = and k — 3. 



2. Stretching interpolation technique 




FIG. 5. The correlation coefficient CC{e) is evaluated at two 
temperatures for k =0 and 3 in the [12.5 /xs — 50 fj,s] range. 
The maximum, obtained for a relative velocity change of £3 = 
1.86 10"'^, is indicated by the vertical arrow. It corresponds 
to an increase of temperature of 0.8° C. Theory is from Eq. [5] 



In the stretching technique, the coda hk{R, R,t) is in- 
terpolated at times t{l~ e) with various relative velocity 
changes, e, in the [ti — ^2] time window. Ek is therefore 
the e that maximizes the cross-correlation coefficient: 



CCkie) = 



I'^hk [t{l-e)]ho[t]dt 



(4) 



An example of correlation coefficient is plotted in 
Fig. [5l If we assume that ho and hk are stationary 
waveforms^ - and are well described by Eqs.fTH^l we have 
a theoretical estimation of CC: 



CCkie) 



A 



p{Lu)sin (ujet2) — sin (tueti) dio 
uje (t2 - ti) 



Bis) 



(5) 



which in the simple case oi ti — and t2 — T simply 
reduces to: 



B{£) is a random process of zero mean and standard 



CCkie) =A^ 



Iau> Pi'^)sinc (w (e - e^) T) duj 



+ Bie) (6) 



with CO the pulsation, Au the bandwidth, piuj) the power 
spectrum density. The constant A depends on the vari- 
ance of G, noted (G^) and the variance of the additional 
fluctuations, noted (/^), 



(7) 



2ti 



TA 



(8) 



The term containing the sine function is represented by 
the crosses in Fig. [51^61 and the fluctuations (y^B^)) 
around this average are in gray. If the amplitude of 
the sine function is much greater than the fluctuations, 
A > 7(^2)", the maximum of the cross-correlation 
coefficient CCk is obtained for e — ek, which provides 
the relative velocity change for the given state k. It is 
interesting to note that the peak of the sine function 
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B. Active experiment: low quality data 

To mimic a practical situation where data includes ad- 
ditional noise (instrumental or electronic) , we add a ran- 
dom (5-correlated noise n{t) of zero mean to the signals 
h{t) in Eq.[T]y2l For simplicity, we neglect the distortion 
/fc in the considered time-window, and assume a station- 
ary noise: 

{nl) = (nl) = (n^) (9) 
then we get a similar expressions as Eq. [THHl with: 



(10) 



and 



FIG. 6. The correlation coefRcient CC{e) is evaluated at two 
temperatures k = 0&3 in the [12.5 /is — 50 /is] range, for var- 
ious signal-to-noise ratio SNR. Crosses indicate experimental 
data. The gray background indicates expected fluctuations 
(Eq. Ill|l around the theory (Eq. llOp . The proper velocity 
change £3 is found in all cases, though a slight difference is 
visible for SNR=1. 



is visible even if the distortion f{t) (or electronic noise 
noise n{t), see next subsection) are strong. In such a 
case, increasing the integration time T or the frequency 
bandwidth Aw can reduce the fluctuations B. This is a 
crucial advantage of the present technique compared to 
the doublet technique in the case of noisy or distorted 



2tt v/("^) +2(^2) (/i2 



TAu 



(11) 



The velocity change is measured again for signal-to- 
noise (SNR) ranging from 1 to 100 (Figs. [71). For a SNR 
of 100 and 10, we find the same results for the stretching 
and for the doublet technique. However, if the SNR is de- 
creased to 2, the velocity variations measured from the 
doublet technique are not accurate at all, while they re- 
main relevant with the stretching technique. This estab- 
lishes the stretching technique as a more stable process- 
ing procedure for noisy records. Note that the connection 
between the fluctuation of the waveforms and the error in 
the estimation of dV/V = e will be subject to a separate 
communication. Nevertheless, this error can be visually 
estimated by the gray area around the theoretical curve 
(crosses) in Fig. El 



C. Advantages and drawbacks of both techniques 

The doublet technique (CSMWT) has been used 
successfully for more than 20 years to efficiently retrieve 
small velocity changes in the mediumii^iii^. This tech- 
nique does not suffer from change in amplitude of the 
waveform, including the coda decay, and the processing 
is very fast. It also manages clock errors in origin time 
without further processing, which is a central issue in 
active and passive field experiments^^. It also allows to 
select a given time window in the data. 

The stretching technique is more recent. It is based 
on a grid-search for e, and is found to be slightly more 
time consuming in terms of computer processing. A no- 
ticeable disadvantage of this latter technique is also that 
it assumes a linear stretching of the waveform, which 
is not valid for media with heterogeneous changes (in- 



cluding the Earth). The main interest of the stretching 
technique versus the doublet one is its stability toward 
fluctuations (noise) in the data, as mentioned in the pre- 
vious section and demonstrated by fig. [71 This provides 
an opportunity to increase the sensitivity of detection of 
weak changes in the earth's crust with seismic w&ybsM^^. 

IV. ACTIVE AND PASSIVE EXPERIMENT: 
MONITORING WITH THE CORRELATION? 

Most previous authors suggested that monitoring 
weak changes in the earth with ambient seismic noise 
correlation is based on the assumption that those corre- 
lations yield the actual GF of the medium. Thus, the 
late part of the correlation is interpreted as the coda of 
the reconstructed impulse response. Is this assumption 
actually necessary to monitor the changing earth with 
good accuracy? We address the question in the present 
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Stretching technique Doublet technique 

FIG. 7. Relative velocity changes evaluated from the stretching (left) and for the doublet (right) technique for various SNR. 
The actual dV/V is satisfiyingly retrieved for any SNR> 1 with the first technique, not with the second. 



section by comparing relative velocity changes measured 
either in the autocorrelations of records from distant 
noise sources (the passive experimental setup), or in 
pulse-echo data (the active experimental setup). 

In the passive experiment, transducers are attached at 
opposite sides of the gel (figure [2]-right). On one side, 16 
sources (S) are acting to mimic a distribution of noise 
sources. Impulse responses are recorded by the 7 re- 
ceiving transducers (R) and consecutively convolved by a 
white noise to mimic acoustic (or seismic) ambient noise. 
Then they are autocorrelated at each receiver. Note that 
the precise knowledge of the noise source position is un- 
necessary: the source position has no effect on the veloc- 
ity change estimation. If the noise sources were uniformly 
distributed and the coda records long enough, these cor- 
relations averaged over sources should result in the GF 
for the medium^. However, the records used in this ex- 
periment are of finite duration, and the auto-correlation 
has not converged to the Green function yet. This can 
be seen visually in figure [SJ where the autocorrelation is 
plotted alongside the (time derivative of the) reference 
Green function obtained in the active experiment. The 
fact that the two signals are uncorrelated is confirmed by 



the low value of the coherence between them (2%). 
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FIG. 8. Comparison of the pulse-echo data h[R, R, t) ob- 
tained in the active experiment and the average autocorrela- 
tion dth{S,R,t) X h{S,R,t). 



The noise signals from each source are emitted at con- 
secutive times. To emulate a signal coming from multiple 
sources at once, the signals recorded from each source i 
are stacked: 



Stability of noise-based correlation monitoring 6 



/ifc(i) = ^Gfc(5„i?,i)®n,(t), (12) 

i 

where the subscript k holds for temperature. The veloc- 
ity variations are then computed using the autocorrela- 
tion of this total signal. 

Ck{t) = hk{t) X hk{t). (13) 

The velocity variations for each execution k of these 
experiments are displayed in figure [9l alongside those 
found with the active experiment. The acquisitions run 
over about 25 minutes, over which the temperature has 
increased by about 0.8°C. Within the errors, the velocity 
variations found with the autocorrelation are the same 
as those found with the reference GF. 

Until now, the analysis was based on the assumption 
that the autocorrelation used in the passive experi- 
ment closely resembles the GF of the medium if there 
are enough sources, and these sources are stable. In 
our experiment, the GF is not reconstructed in the 
auto-correlation. Nevertheless, figure [H demonstrates^^ 
that it is still possible to retrieve correct information 
about small changes in the medium properties with 
the resulting auto-correlation. This is a very promising 
observation that supports the idea that correlation of 
seismic noise will give the opportunity to monitor weak 
changes in the Earth with good reliability even when 
the correlations have not converged to the GF. Indeed, 
in both active and passive experiments, we measure 
the acoustic/seismic signatures of the medium, which 
naturally include its weak variations. 

Nevertheless, in order to achieve a proper comparison 
between our laboratory-scale experiment and seismology, 
we have to take into account another phenomenon. In- 
deed, on the earth, the seismic noise sources location is 
smoothly changing from one week to another. The ques- 
tion addressed in the next section is then: what will hap- 
pen to our monitoring technique when the background 
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FIG. 9. Relative velocity changes evaluated from the active 
(pulse-echo with S=R) and the passive (autocorrelation with 
St^R) setup. 



Noise Sources 




FIG. 10. (a) The seismic noise source structure changes from 
one date to another, (b) Analogous laboratory experiment; 
two sets of sources are chosen at two different dates. 



noise is no longer stable, ie. the source distribution 
changes spatially? 

V. INFLUENCE OF NOISE SOURCE STABILITY 

The change in background noise structure is simulated 
by averaging the autocorrelations for a number of simul- 
taneous uncorrelated sources i. The same is done for a 
slightly different selection of sources a few instants later, 
then we calculate the relative velocity change dV/V be- 
tween the two autocorrelations. A simple picture of this 
is given in Fig. [TOl Imagine for instance that at date 
k =0, sources 1 and 2 are active, and at date k —1 
sources 2 and 3. For three sources, the decorrelation 
of the signals can be analyzed theoretically. For the first 
experiment (fc =0), the record is: 

hQ = G{Si,R,t)®ni{t)^G{S2,Rt)®n2{t) (14) 
and its auto-correlation reads: 

Co = ft-o X ho (15) 

or: 

Co = G(S'i , i?, t) X G{Si ,R,t)® ni [t) x ni{^t) (16) 
-f G(52, R, t) X G{S2,R, t) ® n2{t) X n2{-t){17) 
+ 2G{Si,R,t) X G{S2,R,t)(g>ni{t) x n2(-tOl8) 

Since ni(t) x n2{t) is almost zero, we neglect the third 
term. For simplicity, we shorten the notation: 

G{Si, R, t) X G{S„R, t) ® n,{t) x n,{~t) = AC^{t) (19) 
Which leads to : 

Cn{t) = ACi{t) + AC2{t) (20) 
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FIG. 11. Relative error in the estimation of e = dV/V, versus 
the ratio of unchanged-to-total amount of noise sources (x- 
axis in log. scale). 



Similarly, for the second experiment (fc =1), in which 
the signal is stretched due to the velocity change and the 
sources used are numbered 2 and 3, we get: 

Ck = AC2{t[l - e]) + ACsitil - e]). (21) 

The objective is to find the Sk that maximizes the 
cross correlation coefficient CCk defined in Eq. [H The 
terms ACi and AC3 are decorrelated waveforms that 
play the role of fluctuations (eq. [2]) and will contribute 
to the term B in eq. [5l Assuming that the variance of 
ACk is constant: (ACf) = (ACi) = (ACi), we find 
A — ^ and the standard deviation of B simplifies as 
y/{B^ = V5V2tt/2\/TAlu. From this latter equation, 
we deduce that a proper estimation of dV/V is carried 
out if we process a sufficiently large amount of data: 
VtA^ > VsV2n. This fi gure is valid for two sets of 
sources that have 50% of sources in common. The same 
calculation can be carried out for any ratio of unchanged- 
to-total amount of uncorrelated sources. 

As an example, we report in Fig. [11] the relative 
experimental error in the estimation of dV/V for various 
unchanged-to-total source amount ratio. The velocity 
evolution retrieved through the autocorrelation, as 
shown in figdl is used as a reference. For different ratios 
of spatially unchanged sources, the velocity changes 
are computed again, and the deviation with respect 
to the reference is considered the error. In fig [Tl] the 
relative error estimation is shown for different ratios 
of unchanged-to-total sources. When half the sources 
remain stationary, the relative error is ~ 20 %, meaning 
that we have access to a rough (but relevant and 
interpretable) estimate of the velocity change. With the 
given experimental coda duration and bandwidth, we 
observe that a satisfying estimation of dV/V is obtained 
if 50% of the sources are unchanged. 

In this framework, we conclude more generally that 
the spatial instability of the source distribution is not a 
limitation for noise-based correlation monitoring as long 
as at least part of the noise spatial distribution is stable. 



The smaller the stable area, the harder the dV/V estima- 
tion, and similar to the conclusion of section III, a larger 
the integration time T or bandwidth Alu is required. 



VI. DISCUSSION AND PERSPECTIVES 

In this paper we conducted laboratory experiments 
with ultrasonic waves to monitor weak velocity changes 
in the medium. To that end, we employed and compared 
several procedures that process small phase shift in the 
diffuse coda waves. These phase shifts correspond to 
change of arrival time in the waveforms, which were due 
to small temperature changes in the medium. 

The paper began with an active (pulse-echo) exper- 
iment (section III), in which we compared the doublet 
(or CSMWT) technique to the more recent stretching 
technique. The former is based on Fourier analysis in 
multiple time windows. The latter is based on the in- 
terpolation of the whole waveform and on a grid search 
optimization. The latter was found to require more com- 
puting power, but was also found to be more stable to- 
ward noise in the data. In the second part of the paper 
(section IV), we processed the auto-correlation of noise 
records acquired in the same medium. Active (pulse- 
echo) and passive (auto-correlations) data were processed 
using exactly the same processing procedure: the rela- 
tive velocity changes of the medium were deduced from 
the late arrivals using the stretching technique. Very 
similar results were found in both cases, although auto- 
correlations had not at all converged to the GF. We there- 
fore demonstrated that, contrary to prior belief, pas- 
sive monitoring with ambient noise remains possible even 
when the correlation has not converged to the GF. In 
other words, noise-based monitoring requires weaker as- 
sumptions than noise-based imaging. In the last part of 
the present document (section V), we tested the robust- 
ness of the noise-based monitoring technique in the case 
of unstable distributions of noise sources. We demon- 
strated that, as long as a certain portion of the sources 
is stable, velocity variations can still be retrieved. 

Even though we consider a laboratory experiment in 
this paper, in practice the results can be extended to 
different scales, among which seismology is of particular 
interest. Note that the idea that the coda of the 
correlations contains precious informations has been 
recently demonstrated by studying the Correlation 
of the Coda of the Correlation (C3)2S. In previous 
monitoring studies along the San-Andreas faultline at 
Parkfield, California^, there were some doubts as to 
whether or not the GF was properly reconstructed in 
the correlations. This was due to i) the imperfect source 
distribution in the ocean, ii) to the short time-series 
over which correlations were performed and iii) the low 
quality data for frequencies below 1 Hz. Nevertheless, it 
was still possible to observe small variations in the coda 
of the correlations. From our laboratory experiment, we 
can confirm that these latter changes are actual physical 
observations that can be interpreted by velocity changes 
in the crust. On a somewhat smaller scale, passive 
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monitoring can be applied to seismic prospecting on 
reservoirs. On such reservoirs, the technique could for 
instance be used to follow the effect of fluid flows as oil 
or gas. During production, we indeed expect relative 
changes of velocity in the reservoir of the order of a few 
percent^, which seems to be observable by the method 
presented here. It might also be possible to detect the 
velocity variations caused by stress changes following 
subsidence. This monitoring could either be performed 
with reproducible active sources or with background 
seismic noise. 
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